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“Since ’tis Nature’s law to change, 

Constancy alone is strange” — Earl of Rochester. 

Abstract 

While analyzing second order thermodynamical phase transitions, Lev 
Landau (the famous Russian physicist) introduced a very vital concept, 
the concept of an ’’order parameter”. This not only amalgamated the 
previous fragmentary theoretical understanding of phase transitions (an 
arsenal of mean-field theories) but also it put forward the important theory 
of ” spontaneous symmetry breaking” . Today, order parameter concept is 
a paradigm both in condensed matter physics and in high energy physics, 
and Landau theory is a pinnacle of all mean-field theories. Mean field 
theories are good qualitative descriptors of the phase transition behavior. 
But all mean-field theories (including Landau’s theory) fail at the criti- 
cal point (the problem of large correlation length). The problems with 
large correlation length in quantum many-body systems are the hardest 
problems known in theoretical physics (both in condensed matter and in 
particle physics). It was Ken Wilson’s physical insights and his power- 
ful mathematical skills that opened a way to the solution of such hard 
problems. 

This manuscript is a perspective on these issues. Starting with sim- 
ple examples of phase transitions (like ice/water; diamond/grapliite etc.) 
we address the following important questions: Why does non-analyticity 
(sharp phase transitions) arise when thermodynamical functions (i.e., free 
energies etc) are good analytic functions? How does Landau’s program 
unify all the previous mean-field theories? Why do all the mean-field theo- 
ries fail near the critical point? How does Wilson’s program go beyond all 
the mean-field theories? What is the origin emergence and universality? 
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1 Introduction to Thermodynamical Phase Tran- 
sitions 

Mankind is familiar with various phases of matter and transitions between the 
phases from antiquity. For example, ice-water-steam; ice is solid that exhibit 
characteristic rigidity, water is liquid and it flows like other liquids, and steam is 
a gas, much like other gases. It is only in the 19th century that with the advent 
of the science of thermodynamics the phases of matter and their transformations 
were rationalized within the scientific domain. The phase transition of water 
to steam gave birth to the great discovery of the steam engine by Thomas 
Newcomen in 1712, and it resulted in the subsequent great industrial revolution 
first in Britain and then in Europe. It is beyond doubt that man is able to use, 
for the service of humanity, the great practical potential of the fundamental 
understanding of the materials and their transformations. 

Aim of the the present paper is to give a broad overview of the phase tran- 
sition phenomena from fundamental perspective starting from early works of 
Andrews and van der Waals. The Mean-Field Theories (MFTs) were devel- 
oped to understand thermodynamical phase transitions from atomic/molecular 
perspective, and van der Waals was the pioneer in this [I]. We discuss various 
Mean-Field Theories (MFTs) and their reformulation by Lev Landau in 1937[2j. 
MFTs are good qualitative descriptors of the phase transition behavior. But 
all mean-field theories fail at the critical point including the Landau’s formu- 
lation. At the critical point correlation length (defined precisely in the text) 
diverges and fluctuations dominate over the average behavior. It is this failure 
of the MFTs that triggered further investigations that resulted the formulation 
of renormalization group by Wilson and others in 1970s[3]. The problems with 
large correlation length in quantum many-body systems are the hardest prob- 
lems known in theoretical physics (both in condensed matter and in particle 
physics) and this branch of physics is far from being closed. 

The present paper is a pedagogical introduction to this vast but coherent 
topic. The story of the successes and failures of MFTs and their systematic 


2 


replacement at the critical point by the renormalization group theory is given 
in a clear and pedagogical way. We also address, more specifically, the follow- 
ing questions: Why does non-analyticity (sharp phase transitions) arise when 
thermodynamical functions (i.e. , free energies etc) are good analytic functions? 
How does Landau’s program unify all the previous mean-field theories? Why 
do all the mean-field theories fail near the critical point? How does Wilson’s 
program go beyond all the mean-field theories? What is origin emergence and 
universality? 

The answers of above questions are well known [31 [31- Presentation here 

is sufficiently self-contained to address the above questions. Thus it will be 
useful to the students learning the subject. We also express our standpoint re- 
garding the failure of MFTs, and point out that MFTs faithfully reproduce 
diverging fluctuations at the critical point and divergence of fluctuations can be 
understood most clearly by recognizing the fact that the property of Statisti- 
cal Independence (SI) is violated in systems with long correlation length (when 
various parts of the system become statistically correlated) [6] . And the funda- 
mental thesis of Boltzmann- Gibbs statistical mechanics-fluctuations converge 
as oc jj in sum-function observables [B] -is no longer applicable in systems with 
large correlation length (of the order of the system’s size) making the fluctu- 
ations to dominate over the average behavior. The reason why MFTs predict 
wrong values of critical indices is given in [7], that also clarifies a common mis- 
conception. 

The paper is organized as follows. In the next subsection we give some exam- 
ples of phase transitions. The next section (section 2) is devoted to Ehrenfest’s 
classification and to some mean-field theories. Then in section 3 we present 
Landau’s formulation of MFTs. Section 4 deals with the problems of MFTs. 
After recognizing the difficulties when dealing with problems with long correla- 
tion length we present Kadanoff’s construction and Wilson’s Renormalization 
Group (RG) formulation that rectify the problems of MFTs (section 5). This 
is illustrated by a specific example of </> 4 - theory. Although there are many dif- 
ferent formulations of renormalization group method and has been applied to 
many different problems[3l 31 E) S El IQ, [IT) , our point here is to present the 
essence of the method for pedagogical purposes. In section 6 we also address 
the important question: Why do all the mean-field theories fail near the critical 
point and how does Wilson’s program go beyond MFTs? A short biographical 
sketch of Ken Wilson is given in the appendix. 

1.1 Some examples 

Phase transitions are ubiquitous in Nature! Most common phase transition is 
that of water and ice (Figure 1(a)). We see this when we put liquid water 
in the refrigerator for making ice, a common experience. We know that to 
make ice from water we need to extract some heat for a given mass of water. 
This is called the latent heat and for water it is 334 Joules/gram . Thus to 
make ice from water you have to extract this much amount of heat per gram 
from water (i.e., by cooling it) and to convert back water from ice you have 
to supply that amount (i.e., by heating the ice). You also note that when you 
cool water to convert it into ice, it expands, i.e., the volume of ice for a given 
mass of water is more than the volume of liquid water of that mass, or ice 
floats on water. These facts can be quantitatively described by introducing 
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the thermodynamical free energy [55]. Discontinuity of its first order derivatives 
with respect to thermodynamical parameters leads to above discontinuities. We 
will explain this in the section devoted to Ehrenfest’s classification. 




(a) (b) 

Figure 1: (a) Three phases of water, (b) Phase diagram of a simple substance 
like water. 

Figure 0Jb) depicts a typical phase diagram of a “simple” substance like wa- 
ter for example (only change one has to make is that of a negative slope of the 
liquid-solid phase line for the case of water). Phase diagram depicts the kind of 
phases that are taken by a given substance under a given set of external con- 
ditions. These external conditions are the thermodynamical parameters of the 
system, for example, Pressure (P) and Temperature (T) for a simple substance, 
or, temperature and external magnetic field for a magnetic substance, like, for 
example, a ferromagnetic. 

The phase diagram, as in figure 1(b), is a two dimensional projection on 
the P-T plane of a three dimensional picture, with free energy along an axis 
coming out of the plane of the paper (for example, referring to figure Q2b)) and 
perpendicular to the P-T axis. So in this 3-D picture one has three free energy 
surfaces. A substance is stable in a given region of the phase diagram if its 
free energy surface takes the minimum value. Consider the gaseous regime, in 
this, the free energy surface for the gaseous phase takes the minimum value as 
compared to two other free energy surfaces that of liquid and solid. The phase 
boundary lines shown in figure [T[b) are the lines at which above mentioned 
surfaces cut each other. And at each phase line free energies of two substances 
separated by that line takes the same value. Thus on each phase line both can 
co-exist. The point at which all can co-exist is called the triple point. Another 
important point is called the critical point. It is a point in a phase diagram 
where first order phase transition goes to second order phase transition. The 
discontinuity in the first order derivatives of g (latent heat and change of volume 
per unit mass) goes to zero. Thus in liquid-gas transition the density difference 
(between gas and liquid) becomes zero and it is not possible to distinguish 
gas from liquid at the critical poinlQ. However, the second order derivatives 

1 But density fluctuations dominate over all length scales at the critical point. Density 
fluctuations on the length scale of the wavelength of visible light cause strong scattering of 
light and results in a very spectacular phenomenon called the critical opalescence. Here is 
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of the Gibbs free energy remains discontinuous (see the section on Ehrenfest’s 
classification) . 



(a) (b) 


Figure 2: (a) Allotropes of Carbon: Graphite is soft! while 

Diamond is the hardest material known (Figure source: 
http:/ /mrsec. wisc.edu/Edetc/nanoquest/carbon/index.html). (b) Phase 

diagram of Graphite. 

Diamond and Graphite: Another example of phase transitions is that of di- 
amond and Graphite (figure ([%])). These are two allotropes of Carbon. Pencils 
write because layers of graphite easily slip past each other and can coat an- 
other object, such as paper. Diamonds look very very different from pencils! 
They are the hardest mineral known!! This is an example of the structural 
phase transition in which Carbon atoms rearrange themselves. Diamonds and 
Graphite are naturally occurring substances. Diamonds formed from Carbon 
deep under earth (~ 200 km) by geologic processes under larger pressure and 
temperature environments (as compared to ambient conditions). One can even 
transform Graphite into Diamonds-in a Lab! This is possible by applying a 
very large pressure, 100,000 times that of normal atmospheric pressure, with 
graphite sample kept at 2000 K (figure [21 (b)). Thus formed artificial diamonds 
are smaller in size as compared to diamonds created by geological processes. 
Nature is a smarter creator! One would like to say! 

Superfluidity and Superconductivity: Phase transitions can also lead to 
quantum behavior at macro-scale. Superconductivity and superfluidity are two 
prime examples [ 131 . One interesting example is that of Helium. Helium exists 
in nature in two stable isotopes. One common isotope is 2 He 4 . The other rare 
isotope 2 He 3 is one atom for a million atoms of 2 He 4 in earth’s atmosphere. 
Consider first the phase diagram of common isotope 2 He 4 (Figure EKa))- It liq- 
uefies at 4.2 K at the ambient pressure (Heike Kamerlingh Onnes was the first 
to liquefy helium. He liquefied it in 1908 and later discovered superconductivity 

a quick experiment 1121 . Take equal quantities of cyclohexane and acetic anhydride in a test 
tube. Two liquids do not mix at room temperature, but on heating they mix. If one let 
the mixture to cool, near 53 C the mixture turns milky white — -A spectacular phenomenon 
called critical opalescence. 
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Figure 3: Phase transitions in Helium, (a) Common Isotope 2 He 4 , (b) The rare 
isotope 2 He 3 . 


in Mercury in 1911). When cooled further ~ 2.17 K it transforms to a new 
form of liquid (a new thermodynamical phase called the superfluid phase). In 
this phase liquid Helium flows without friction (zero viscosity) and becomes a 
“super- heat-conductor” . This superfluid phase of 2 He 4 is also called the Hell 
phase (the normal one is called Hel). Heat capacity shows an infinite disconti- 
nuity at 2.17 K (figure (fT]l). There is no latent heat involvec0. The shape of 
heat capacity curve resembles with that of a Greek letter lambda (A). Thus it 
is also called the lambda-transition in liquid Helium literature. Liquid Helium 
can be made solid by applying a pressure of 25 atm for temperature less than 2 
K. Higher temperatures require higher pressure for solidification (figure [3] (a)). 
Interestingly it does not have a triple point! The rare isitope 2 He 2 3 behaves 



Figure 4: The lambda-transition in common Helium, 
exotically (figure EDb)). If you start from a very low temperature, say 0.0001 K, 

2 But it is not the second order phase transition! see Ehrenfest’s classification in the next 

section. 
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Figure 5: Superconducting phase transition, (a) Heat capacity show a finite dis- 
continuity at T c , (b) Superconducting substances show Meissner effect (Image: 
public domain/ Wikimedia Commons.). 


with a sample of 2 He 3 at 30 atm pressure, you first have liquid Helium. On 
reaching 0.18 K it solidifies! It remains solid upto 0.49 K and again turns liquid 
for higher temperatures. This behavior of the rare isotope of Helium is quite 
unusual (as common solids remain solid with reducing temperature!). 

The superfluidity of the normal isotope 2 He 4 is due to the Bose-Einstein 
condensation below the condensation temperature ~ 2.17 K. As 2 He 4 atoms 
are bosons they condense into the ground state with order parameter (the con- 
densate wavefunction) as pe^. Where p is the superfluid density and <p is phase 
of the wave function. All the observed properties can be understood using two- 
fluid picture (see for details [H]). The reason why Helium remains liquid down 
to the lowest temperatures achievable at normal pressure is the large amplitudes 
of the zero-point quantum oscillations that disrupt the crystalline order. 

The case of the rare isotope is much more interesting. As 2 He 3 atoms are 
fermions they cannot show BEC directly and thus superfuildity. But experimen- 
tally superfluidity has been seen. The reason for superfluidity in the rare iso- 
tope is the formation of Cooper pairs between Helium atoms much like the BCS 
(Bardeen-Cooper-Schrieffer) theory of superconductivity [15], but with impor- 
tant differences: As there is no crystalline lattice, phonons are not operative in 
forming Cooper pairs (in contrast to simple metals explainable by BCS theory). 
The origin of Cooper pairing in the rare isotope is magnetic and Coopers pairs 
resides in higher angular momentum states with internal degrees-of-freedom (see 
Anderson- Morel [T()| ) . 

Superconductivity is an another example that show macroscopic quantum 
effects. First discovered by Heike Kamerlingh Onnes in 1911, it took about 
fifty years of theoretical efforts and finally the microscopic understanding came 
with BCS theory of superconductivity in 1957[T5]. Loosely speaking the basis of 
superconductivity is the Bose-Einstein condensation of Cooper pairs formed due 
to mutual attractive interaction mediated by the exchange of lattice phonons. 
The important difference between rare isotope of Helium and superconductivity 
in simple metals is that the size of a Cooper pair in metals is much greater than 
mean particle separation (~ 1000 times the lattice constant) while in Helium 
it is like forming a bosonic molecule by two fermionic atoms (the size of the 
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Cooper pair is mean particle separation). BCS theory quantitatively account 
for many observed properties of simple superconductors (excluding Cuprate 
superconductors [T3> 'T70. 

Magnetism of Iron: Iron loses its magnetic properties above 770 C (the 
Curie temperature) . What happens is that due to thermal agitation the ordered 
arrangement of electronic spins becomes unstable. The origin of magnetism in 
Iron is due to a dual role played by d-band electrons [X51 IIP] . They participate in 
metallic bounding as well as act like localized spins. Iron drives its magnetism 
from electron-electron repulsive Coloumb interaction due to which they avoid 
each other in real space and that is effectively executed by keeping spins in the 
same direction. This a rough picture (see for more accurate and quantitative 
analysis 1X5]). 

Cosmological Phase Transition: The origin of SU3, Isospin, and Parity! It 
has been postulated (after the discovery of Hubble expansion) that when our 
universe was very young (only sub-pico seconds of age) after the Big-Bang, 
our universe underwent a phase transition from that high energy scale-called 
the Cosmological phase transition. Below that high energy scale the electroweak 
symmetry is broken. This the phase transition of vacuum itself in which vacuum 
expectation value of a scalar field takes non-zero value [2D]. This is intimately 
connected to the Higgs mechanism [2X1. 

The above are some examples of phase transitions. Obviously, there are 
many more examples. The above are sufficient for our purpose of illustrating 
the mean-field theories and renormalization group. Before we discuss the MFTs 
we first review the work of Ehrenfest on the classifications of phase transitions. 


2 Phase transitions without Landau’s paradigm 


2.1 Ehrenfest’s classification of phase transitions 

Paul Ehrenfest (1880-1933 jf| classified various phase transitions in 1933 accord- 
ing to the behavior of the derivatives of the Gibbs free energy (unfortunately, 
later that year he killed himself as he was suffering from severe depression). 
According to this scheme, the order of a phase transi- 
tion is given by the order of the discontinuous deriva- 
tive of Gibbs free energy. On the line separating 
the two phases (say, in the liquid-gas transition) the 
Gibbs free energies per unit mass g(T , P) of the two 
phases take the same value. On the left of the line, g 
is less for the liquid, and on the right of it it is less 
for the gas. By definition [221: g(T,P) = U — Ts + Pv 
or dg = — sdT + vdP with 



Entropy /mass = s = — 


Figure 6: Ehrenfest and 
his son with Einstein 
1920. 


dg 


Volume/mass = i> = — — i t 
oP 

First order phase transitions are those in which first order derivatives of the 
Gibbs’ Free energy are discontinuous: i.e. , if the above derivatives are discon- 


3 He was a student of Ludwig Boltzmann. 


tinuous. If we cross the liquid-solid line vertically in up direction (at constant 
temperature, see figure [T| (b)) clearly: 


dgi , ^ dg s 
~dP T > ~dP' T 


As the system goes from liquid stability to solid stability (an experimental fact!). 
This implies that: 


vi = 


dgi 

dP 


T > V s 


dg s 

dP 


In other words, for a given mass, volume of the liquid phase is more than that 
of the solid phase, or, substances expand on melting^ 

If we cross the liquid-gas line horizontally (ref. figure [Ijb)) from left to right 
(at constant pressure) clearly: 


dg s 

dT 


p > 


dgi | 
dT 


As system goes from solid stability to liquid stability (again an experimental 
fact!). Above implies that 

Ss > - Si 


or 


s s < si , 




qi - g s 

T 


> 0. 


Where qi — q s = l is called the latent heat of solid-to-liquid transition, which 
is positive! as it should be! Thus heat must be absorbed from solid to liquid 
transition. Examples of first order are: Solid-liquid- vapor, transitions in simple 
fluids, allotropic transformations (like diamond graphite) etc. 

For the second order phase transitions, first order derivative of Gibbs free 
energy, i.e., v. s, are continuous, but the second order derivatives, i.e., heat 
capacity c, and compressibility n are discontinuous. Thus no latent heat and 
volume change occur during a second order phase transition. But in a second 
order transition, heat capacity and compressibility show finite discontinuity. An 
ideal second order transition is the superconducting transition in zero magnetic 
field [22]. 

Phase transition in liquid Helium is not second order (as the discontinuity 
in heat capacity is infinite). This is a special transition called the A— transition. 

Above is called the Ehrenfest’s classification scheme. One can define higher 
order phase transitions along these lines (third and higher order derivatives of 
the Gibbs free energy), but those are seldom observed in practice. Ehrenfest’s 
scheme is based on the fundamental thermodynamical principle: the phase of 
minimum free energy is the stable one |22| . From this it derives other properties 
of phase transitions like latent heat, volume changes etc. Importantly, it does 
not provide more fundamental statistical mechanical formulation of the phase 
transition behavior. We will see in the next section that MFTs are the first 
steps in these directions. 

We will also see that with the Landau’s paradigm above scheme of Ehren- 
fests’ can be nicely re-interpreted by using the order parameter concept. 


4 Note that the case of water is opposite (ice is lighter) as the ice-water (solid-liquid) line 
has negative slop. 
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2.2 The Mean-Field Theories(MFTs) 

The above is the macroscopic way of understanding/classifying the phase tran- 
sitions without any microscipic (atomic/molecular) point of view. Attempts 
to understand the laws of heat from the microscipic (atomic/molecular) per- 
spective resulted in the science of statistical mechanics and Ludwig Boltzamnn 
(1844-1906) was the pioneer [23] . Similar to this philosophy, to understand phase 
transition phenomena from atomic/moleculaer perspective, early attempts were 
made by van der Waalfl 

Thus the credit of first mean- field theory goes to van der Waals’jTj. Van der 
Waals was interested in the notion of continuity of matter, namely, how a gaseous 
substance, on cooling, transforms into a liquid while chemically, at molecular 
level things remain intacfl His main concern was the ideal gas equation and 
the problems of the Boyle’s law. Here is a rough sketch of his argument. 

We know that the ideal gas equation 

PV = Nk B T = nRT 

shows no phase transition for obvious reasons (due to the neglect of inter- 
molecular forces). In contrast, real gases at low temperature and high pressure 
deviate markedly from the ideal gas behavior. To see this, let be the molar 
volume (volume per mole) of an ideal gas: 

pV° = RT. 

Define the compression factor Z at a temperature T for a real gas ( PVm = RT): 

ry pVm 

In figure EJa) compression factors of both the gases are plotted as a function 
of pressure. In an ideal gas, as there is no interaction between the molecules, 
the compression factor remains unity with pressure while for a real gas it first 
decreases and then at higher pressure it starts to rise again (figure EK a)). This 
can be easily correlated with the inter- molecular potential plotted in figureEKb). 
At high pressure gas is compressed and molecules are relatively near to each 
other and thus exert repulsive forces on each other as the potential is large 
positive (figure Elb))- This makes the compression factor greater than that of 
an ideal gas. In other words, for the same number of moles of a real gas and 
an ideal gas in the high pressure limit (when the mean separation between the 
molecules is less than r c ), same pressure will be exerted only if the real gas 
occupies more volume thus making Z > 1. In the opposite limit of low pressure 
(mean inter molecular separation greater than r c ), molecules exert attractive 
forces and thus smaller molar volume and smaller compression factor for the real 
gas. To account for deviations of the ideal gas law and to understand the phase 

5 The thought process of understanding macro-phenomema from microscopic perspective 
is very old one, and is related to one of the very basic aspect of human psychology, namely 
the habit of understanding phenomena around from the basic “building blocks”. Examples 
go back to BC era, for example, Democritus was the first to postulate atoms and atomic 
constitution of matter. Other example is that of statistical mechanics (thermodynamics from 
atomic/molecular perspective) and so on. And this unique habit of mankind resulted in a 
body of organized knowledge that we call science. Author do not want to go into the debates 
of reductionism and emergence philosophies [24] 

6 There is no chemical reaction when water converts from liquid to steam and vice versa. 
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Figure 7: (a) Behavior of the compression factor for an ideal gas as compared 
with a real gas. (b) Inter molecular force as a function of separation in a real 
gas. 



Figure 8: Comparing isotherms obtained from van der Waals’ equation with 
those of experiment: curve 1 for T > T c ; curve 2 at T = T C \ curve 3 at T < T c . 
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transition of gas-to-liquid, van de Waals postulated a modified equation of state, 
now called the van der Waals equation of state [I]. This is a remarkable piece of 
physical intuition and deduction of an equation of state from average properties 
of molecular interactions without going into formal mathematical development. 
This is a physical tour de force and he was awarded with Nobel prize (1910) for 
this. 

He argued “....as you are aware the two factors 
which I specified as reasons why a non-dilute aggre- 
gate of moving particles fails to comply with Boyle’s 
law are firstly the attraction between the particles, 
secondly their proper volume...” -From Nobel prize 
lecture— JOHANNES D. VAN DER WAALS. 

Thus replacing the actual volume and pressure in 
the ideal gas equation by an effective ones: 

Nk B T 

Peff = 


V r 


P = 


eff 

Nk B T ,N. 2 

v^m- a{ v ) 



Figure 9: van der Waals 

results a mean-field theory. This takes “an average 1923)). 

effect” of molecular interactions, and this type of theory is called the mean-held 
theory. As the real molecules repel at short distance, the effective volume is less 
V e ff = V — Nb (where N is the number of molecules and b is an effective volume 
of a single molecule). Similarly effective pressure P e // is less due to long range 
attarctive interaction between the molecules. If we compare the isotherms of 
a real gas with those obtained by van der Waals equation we get a qualitative 
agreement. Van der Waals isotherms look qualitatively the same as those in the 
experimental one (figure ©), but this has two basic defects: 

1. Unphysical van der Waals loops. 

2. Different critical point behaviour. 

The unphysical loops means that if one increases volume pressure increases (for 
the curve with positive slop in figure © in the coexistence regime). This is 
unphysical and can be rectified with Maxwell constructiorQ- 

The critical pressure, temperature, and volume can be 
calculated from: 



dp 

dV 


d 2 p 

dV 2 


= 0 . 




And near the critical point \T C — T\ « T c and from the 
above derivatives, it can be easily shown that: 

Pliquid Pgas OC (T c T)^. 

From van der Waals’ theory /3 comes out to be 1/2 (the 
mean- field value). But, in actual experiments on liquid- 
gas systems, /3 = 1/3. This was first experimentally de- 
duced by T. Andrews in 1869! This very clearly showed 


Figure 10: Thomas 
Andrews (1813- 
1885) 

the contradiction of the mean- field theory with experiment. 

' In Maxwell’s construction one replaces the curve in the co-existence regime by a straight 
horizontal line such that the area under it is equal the area above it (ref. figure [8] (a)) [25]. 
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We will see that this discrepancy is not specific to van der Waals theory but 
a basic defect of all such theories called the MFTs at the critical point. It took 
about 100 years for the resolution of this discrepancy!! This happened in the 
hands of K. G. Wilson in early 70s. The culprit is non-negligible fluctuations 
(density fluctuations in this liquid-gas case) at all length scales that invalidates 
the application of the mean-field theory itself [7]. 

The second example of MFT we consider is that of the famous Ising model — 
a caricature for spontaneous ferromagnetism (like in Iron) . A simple two dimen- 




Figure 11: (a) 2-D Ising model (b) Phase diagram of the model. 

sional Ising model is depicted in figure HU a) in which spins are represented by 
arrows and these interact with their nearest neighbors. Each spin can have two 
states, either up, or down. This short range interacting model shows the phe- 
nomenon of spontaneous magnetization. Consider that the model is immersed 
in an external magnetic field H. Write the Hamiltonian in dimensionless form: 

kwf = ~ K Y 

** <nn> i 

Here I\ = and h = There are two competing tendencies: 

1. Lining- up of spins due to mutual spin-spin interaction ( K ) and the exter- 
nal field [h). 

2. Disruption due to thermal agitation. 

At low temperatures first tendency wins over and at high temperature the second 
one wins over. To calculate mean magnetization (a measure of the lining-up of 
the spins), one computes the Free Energy in the standard way: 


F = — fc^Tln I e ~H(c)/k B T 

\ [all configurations : c] 

Here comes the big question: Why do statistical mechanical systems show phase 
transitions (non- analytical behavior like spontaneous magnetization, see figure 
llll bk the graph near T = T c ) when the system’s free energy is a smooth 
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function of temperature and of system’s parameters (an analytical function)? 
Sum of smooth exponentials is a smooth function! 

The origin of non-analyticity (like spontaneous magnetization) is in the ther- 
modynamical limit (infinite system size!): Larger the system size sharper is the 
transition: As real phase transitions involve thermodynamical systems having 
very large number of particles of the order of Avogadro number, transitions 
appear infinitely sharp to human perception! This has been dubbed as “ex- 
tended singularity theorem” by Leo Kadanoff or better one can say “extended 
singularity conjecture” [2] • As the proof exists only in one special case PHI. 

Mean-field theory for the Ising model is done in usual wav [41. 127: l28|: 

1. First, consider that just one spin is immersed in magnetic field (that can 
take two values: +1 (up-spin, say), and —1 (down-spin)). Thermal expec- 
tation value of the spin is: 

W> = -v = tanh ( ft )- 

2^a=±l e 

2. Second, when our test spin is interacting with many neighboring spins, all 
immersed in magnetic field: 

( a ) = ^ g=±1 lh effa = ta nh (h + I< ^ a t ) ~ tanh (h eff ). 

1 <nn> 


h e f f = h + effective field due to other spins = h + zK (a) . 

Here the exact field K 'f2 <nn> <Ji due to near-neighbor spins is replaced 
by the effective field ~ zK(a) (z is called the coordination number, = 4 
in 2D). This constitutes the mean- field approximation. 

Thus in MFT: (a) = tanh(h + zK(a)). To see what values (a) can take, 
one must solve this self-consistently. As K oc ^ it turns out that zK = 
(where T c is the critical temperature). In zero external field h = 0 the equation 
(a) = tanh(^(a)) has the following solutions (figure [12]): 

f 0 if T>T C 

(a) = l 0, +1, — 1 if T <T C 
y degenerate if T = T C 

The region T <T C is of particular importance. Near the critical point (ex) — 0, 
tanh(a;) can be expanded for small x as ~ x — ^ x 3 + .... And it can be easily 
shown that near the critical point: 

(a) oc (T c — T)P. 

With /3 = f (the mean field value). Compare this with van der Waals theory! 
(THE VERY SAME BETA VALUE!) 

Experiments on a wide variety of magnetic systems show that ff ~ | . Again, 
we see the contradiction between mean-field theory and experiments. So, some- 
thing is wrong with mean- field theory near the critical point! But, as mentioned 
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Figure 12: Solutions of the zero field mean-field equation. 


before, solution comes much later. However, the great merit of MFTs is that 
they qualitatively explain the emergence of phases below T c (refer to the middle 
graph in figure m- One can show that the free energy for (a) = 0 phase is 
greater than the other phases with non-zero (<r), and the phase with minimum 
free energy is the phase that is realized in a thermodynamical system. Thus one 
has non-zero (a) for temperature less than the critical temperature (figure HID . 


3 Landau’s paradigm: symmetry breaking and 
the order parameter concept 

In a landmark paper[2] of 1937 Lev Landau formu- 
lated the phase transition problems from the broken 
symmetry perspective. All the above phase transi- 
tion theories can be written in one mathematical 
language near the critical point. He noticed that 
most of the phase transitions involves breaking of 
discrete or continuous symmetry and introduced a 
very important concept to capture the nature and 
extent of this symmetry breaking: the order param- 
eter concept. One can appreciate this rather formal 
language related to symmetry breaking by consider- 
ing some concrete examples. Before we do that, let 
us first consider two important theorems related to 
symmetry and its breaking. Emy Noether in 1918 
put forward an important theorem: The Noether 

theorem (1918): 

“Every continuous symmetry of a system can be associated with a conserved 
quantity.” This theorem is the cornstone of classical mechanics. There are 
numerous example, consider for example; if, gravitational potential along x- 
direction is constant, then linear momentum along that direction is conserved 
and so on. Another important theorem, for the present context, is the Goldstone 
theorem (1961): “The spontaneous breaking of a continuous symmetry can be 
associated with a spinless and massless particle-the Nambu-Goldstone boson 
(or simply the Goldstone Boson)” . To illustrate this we consider few examples: 
Example No. 1: 

1. liquid-to-Crystal transition: In the liquid state atoms/molecules are un- 
der the effect of mutual potential (like van der Waals potential) which 



Figure 13: Lev. Landau 
(1908-1968) 
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Figure 14: Emmy Noether (1882-1935) and Jeffrey Goldstone (born 1933). 


does not break the translational and rotational symmetry of the space 
(molecules/atoms are free to wander around and rotate). 

2. In the crystalline state, a periodic potential emerges, which break the 
continuous symmetry under translations (however, this leads to a discrete 
translational symmetry due to the periodicity of the lattice: the crystal 
momentum) JT5) . 

3. Since in the above example continuous translational symmetry is broken, 
momentum does not remain a conserved quantity. And this violation of 
the conservation of momentum in the crystalline state leads a Goldstone 
boson called the phonon (lattice vibrations). 

Example no. 2: Ferromagnetic transition and broken rotational symmetry: In 
a ferromagnetic system above the transition temperature, spins are randomly 
oriented in all possible directions. But below the transition temperature, this 
rotational symmetry is broken and spins align in the same direction. The Gold- 
stone boson in the ferromagnetic case is the spin- wave excitations (or, collective 
excitation of spin waves). 

One can consider many more examples. Thus each phase transition is a 
manifestation of some kind of broken symmetry: In ferromagnetic rotational 
symmetry is broken; in crystal formation translational symmetry is broken; in 
superconductor or superfluid electromagnetic gauge symmetry is broken.... On 
this basis, Landau introduced the very important concept of an order parameter. 
For example, in ferromagnetic transition the order parameter is an average mag- 
netization per unit volume (which is zero for T > T c and non-zero for T < T c ); 
in liquid-gas transition it is the density difference; in superconductor it is the 
expectation value of the Cooper field < t/'J-V’-o- >1 and hr electroweak symmetry 
breaking it is expectation value of Higgs field < <f> >.... Order parameter is zero 
in symmetrical phase and is non-zero in the unsymmetrical phase. It is a macro- 
scopic variable, much like a hydrodynanrical variable. It has discontinuous jump 
for first order transitions as the order parameters are the first order derivatives 
of the free energy (for example, M = — J^). Thus one can connect Landau’s 
scheme with that of Ehrenfests’ and one can interpret Ehrenfests’ scheme in 
terms of continuity or discontinuity of order parameters. For the second order 
transitions order parameter evolve from zero in symmetrical phase to non-zero 
in the unsymmetrical phase in continuous manner (thus 2nd order transitions 
are also called continuous transitions). Landau’s great contribution in unifying 
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the theoretical understanding of wide variety of phase transitions lies in his se- 
ries expansion of free energy as a function of order parameter near the critical 
point where the order parameter is small. 

G(M, T) = G 0 (T) + r(T)M 2 + u(T)M 4 + ... 

Where M is some order parameter (like magnetization in ferromagnetic case). 
The fundamental assumption that goes into this is the analyticity of free energy 
as a function of order parameter (magnetization here) and temperature. This 
fundamental assumption gets its justification a-posterior! Most importantly, 
this is in accord with ’’extended singularity conjecture” in which free energy 
per unit volume remains finite in the thermodynamical limit ( N — > oo, V — > 
oo, y = finite). The sharpness of phase transitions (apparent non-analyticity) 
is something to do with very large number, of the order of Avogedro number, 
of particles present in the system. 

Landau argued, according to the basic principle of thermodynamics, that 
at a given temperature, the phase with minimum free energy is the stable one: 
= 0 leads to: 

r(T)M + 2u(T)M 3 = 0. 

Which implies that either M = 0 or M = ± ^j— 

One can relate these to the following two possible cases: 

1. Case A: If r(T) > 0: one has only one minimum at M = 0 ( u(T ) has to 
be positive as free energy is bounded from below) . M = 0 in this case and 
it corresponds to symmetrical phase (T > T c ). 

2. Case B: If r(T) < 0: one has two minima at M = ± yf — and this 
corresponds to the unsymmetrical phase (T < T c ) (see figure (fl5l) l with 
non-zero value of the order parameter. 

From above, one can deduce the critical exponent /? as follows: 

1. As r(T) is positive for T > T c and negative for T < T c , it must pass 
through zero at T = T c . 

2. So near the critical temperature r(T ) can be expanded as: r(t ) = a(T — 
T c ) + 0((T — T c ) 2 ). 

3. Thus magnetization in T < T c phase can be written as: M = ±*J ^ 
0 T-T c y . Here fi = \. 

Again /3 takes the same value 1/2, as in van der Waals theory and in the mean- 
field theory of the Ising model. The general disagreement of MFTs with ex- 
periment is clear. We will see in the next section that systems near criticality 
violates an important property of statistical independence and the correlation 
length diverges at the critical point. Before we embark upon the question of the 
behaviour of correlation length and renormalization group we first would like 
to consider the question of the critical point in the next section in the light of 
some exact solutions. 


17 


4 Problems of the Mean-Field Theories and Lan- 
dau’s paradigm 


In the previous section we have seen the contradiction 
of MFT with experiment near a critical point. The 
P theory = \ and /3 experiment = §• This was apparent 
even in the first decade of 20th century. Still, MFTs were 
developed (like Curi- Weiss, Ising model, Landau’s gener- 
alization) and trusted. Nobody attempted to resolve the 
problems of MFTs at the critical point [2]. Then in 1944, 
with the Onsager’s — a mathematical tour de force — exact 
solution of the 2-D Ising model, it became very clear that 
something is very much wrong with MFTs near critical- 
lity: 

1. C. N. Yang with an extension of Onsager’s solution 
found that /3 is 1/8 different from mean- field result 

P = 1 / 2 - 

2. Onsager’s solution showed that heat capacity di- 
verges as a logarithm of T — T c whereas MFT shows 
a finite jump in heat capacity at T c . 



Figure 15: Behav- 
ior of free energy 
in Landau’s formu- 
lation. 


3. The correlation length £ diverges as -^=== in MFT 

whereas it diverges as T ? T in the Onsager’s solu- 
tion. Thus, v, the critical index of the correlation length, is one half in 
MFT whereas it is equal to one in the Onsager’s solution! 


Here comes the 2nd big question: why do MFTs fail to predict the correct 
values of critical indices? The answer lies in the fact that mean-field equation 
(a) = tanh(h + -p-(cr}) fails to take into account the longer length scale fluctu- 
ations (as compared to inter-particle and lattice constant). See for details)?). 

As fluctuations at all length scales (from lattice constant to the correlation 
length) become important near a critical point, a proper theory in this regime 
must address this issue. The theory is called the renormalization group (next 
section). 


5 Going beyond Landau’s paradigm— The renor- 
malization (semi)group (RG) 

The renormalization theory in general deals with problems of long correlation 
length. The presence of long correlation length leads to the breakdown of an 
important property of statistical mechanical systems called Statistical Indepen- 
dence (SI). Thus the results of standard Gibbs-Boltzmann statistical mechanics 
become unreliable. Near the critical point correlation length diverges as can be 
easily seen: 

Consider that the order parameter is space dependent M(x) and vary very 
slowly as compared to atomic dimensions: One writes Ginzburg-Landau free 
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energy functional [3]: 


F = 


d 3 x{(VM(x)) 2 + r(T)M( xf + u(T)M(a;) 4 - B{x)M(x)}. 


Correlation length £ is defined as the distance over which the effect of a test 
spin extends out to other nearby spins. To calculate £ consider B(x) = B 0 S 3 (x) 
(magnetic field in a narrow region of space). Neglect the smaller 4th order term 
and by minimizing F with a variational calculation, one readily obtains: 

— V 2 M ( x ) + r(T)M (x) = B 0 6 3 (x). 


r(T)|*| 


. From 


With Fourier transforms one gets an immediate solution^/ (x) r 
which one defines a correlation length £ = oc ^ T _ 3 . Thus correlation 

length diverges at the critical point with critical index v = 1/2. 

Due to this diverging correlation length various parts of the system respond 
in a correlated way. This makes the system non-self-averaging! The funda- 
mental hypothesis (actually an important property of ordinary statistical me- 
chanical systems with short range interactions) of statistical independence as 
highlighted by Landau[6] cannot be applied and fluctuations in the additive ob- 
servables (sum-functions) do not obey -^= law(B]. In other words, fluctuations 


could be significant as compared to the average behavior. A comparison of phys- 
ical problems with short and long correlation length is presented in table 0 
The problems with large correlation length are among the hardest problems in 


Table 1: Physical problems with short and long correlation length 


Physical Problems 

Problems with short correlation 
length (£ << L) 

Problems with large correlation 
length (£ ~ L ) 

Various parts of the system with size 
< £ behaves statistically indepen- 
dently 

As correlations extend throughout 
the system, various parts of it are 
statistically correlated! 

Statistical independence (SI) leads 
to product distribution functions: 
Pl2 = PlP2 

P 12 ^ P 1 P 2 

Entropy and Energy of sub-parts be- 
come additive (weakly interacting 
sub-parts of size £) 

non-additivity of entropy and en- 
ergy of sub-parts 

Standard statistical mechanics can 
be applied 

Results of standard statistical me- 
chanics are not reliable! 

Fluctuations in observable converge 
1 

~ vW 

As SI breaks down, fluctuations in 
observable diverge! 

Examples: systems away from crit- 
icallity: MFTs; Hartree-Fock; Clus- 
ter and virial exp... 

Thermodynamical systems near 
criticallity; QCD; Knodo problem; 
cuprate superconductivity etc. 


theoretical physics. For example, thermodynamical systems near criticallity; 
Kondo problem; high temperature superconductivity; relativistic quantum field 
theory; strong interactions and QCD etc. And most importantly, systems with 
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Figure 16: Kadanoff’s Block spin renormalization idea 


large correlation length show a certain kind of cooperative behavior in which the 
details of the microscopic interactions play secondary role! Before we explore 
the renormalization group as formulated by Wilson, we first review the essence 
of a key work due to Kadanoff[55] 0] on which Wilson highly relied. 

5.1 Kadanoff’s great insight 

1. Consider for concreteness the 2-D Ising model. Away from critical point 
(T = T c ) the correlation length is small (£ oc y ), fluctuations on a 
smaller length scale are important and one can reliably use the mean field 
theory. 

2. At T = T c the correlation diverges, as shown before. A reliable picture 
can only be obtained if one take all the spins into account in an area ~ £ 2 
(A HOPELESS TASK !). 

3. Kadanoff approached the problem with a brilliant idea: map the problem 
of large correlation length to a problem with small correlation length!!!! 

4. This is possible by exploiting the fact that on smaller lengths scales l << £ 
spins are well correlated! 

To execute the point number 3 above, consider the original lattice (figure fl6ll 
with lattice constant = a in which we mentally construct blocks, each containing 
four spins (dashed squares in the figure). We then assign an effective spin to 
each block-tlie block spin. The effective spin is an average of the original spins 
(one can use simple majority rule, if majority is up the effective spin or block is 
also up). The ease with which we can define block spins is granted by the point 
4 above. If all the spins are up in a given block, then our effective spin is four 
times the original spin. The resulting lattice of this effective block spins will 
have lattice constant a± = 2 a. The block spins or effective spins will interact 
with each other with some modified interaction that one has to calculate from 
the detailed model. The correlation length scales as follows: If the £ is 1000 a 
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(in the original lattice), then it will be 500 ai (in the scaled lattice). 

Thus the problem with £ = 1000 is reduced to the problem with £ = 500! 
Repetition of this procedure ultimately leads to £ ~ 1 (for a finite system). But 
for a system in thermodynamical limit (system of infinite extent) one can exploit 
the following fundamental fact: “£ diverges at the critical point” This funda- 
mental fact in an infinite system leads to a very curious behavior of system’s 
couplings. They tend towards a fixed point as explained below: 

Let the Hamiltonian of un-scaled system be H = Y^ n » Ko s nS n +i- Here 
Ko is near neighbor interaction in the un-scaled lattice. After first scaling, let 
K 0 — >• and s n Sn' 1 (i.e., block spins Sn^ interact with effective coupling 

K f ^ ) . Thus after first scaling: 

K O s nSn+i K l ffs n )s n+i- 

n,i n,i 


After to repetitions: 

Tr(jn) _ \ ' rs-eff (m) (m) 

11 ef fective / > Iv m ^ n ^ n+i ' 

n,i 

The correlation length behaves as follows: 

After first scaling = |£(iLo) (If the £ is 1000 a (in the origi- 

nal lattice), then it will be 500 ai (in the scaled lattice)). After ?n-scalings 
: = ^Z(Ko). Now the fundamental assumption that goes in 

Kadanoff treatment is this: The scaling does not change the physics of the 

problem, namely, diverges as — >. K c as the original correla- 

tion length £(ATo) does at K c ( K c = k ^ T ). If one knows the functional form 

Kfl f = f{K^£f) for one transformation, the same functional form will be true 
for all scalings! — precisely at the critical point I\q = K c (where correlation 
length is infinite). This scale invariance f figure HTT aB leads to m independent 
relation (AT THE CRITICAL PONIT): 

K c = f(K c .) 

The self consistent solution to this equation leads to the critical point solution 
(or fixed point solution) that gives T c (figure fT7l bll. Thus from the Kadanoff 
picture T c can be calculated if one knows the function f(K). More importantly, 
the critical index v can be expressed in terms of the derivative of f(K ) wrt K 
at K = K c : Suppose K is near to I \ c , then f(K) = f(K c ) + (K — K c )\. Here 
A = As K c = f{K c ), implies f(K) — K c = X(K — K c ). Now suppose 
correlation length behaves as £(A') ~ (K — K c )~ v , then one must have 

£(/(*0) _( f(K)-K c y v 

m v k-k c ) 

From the Taylor expansion above, it immediately follows that \ = (A) - " as the 
scaling factor is 2. Thus if A is known then v can be determined (see for other 
complications [3] ) . 

The problem of the Kadanoff ’s procedure is that it does not give the central 
method to determine the important function f(K). He used the above argument 
to derive scaling relations [33]. Important point worth noting is that Kadanoff 
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Figure 17: (a) Scale invariance: same kind of structure re-appears at all length 
scales when correlation length diverges (if one magnifies any part of the system, 
then it is not possible to notice the difference between the magnified part and 
the original system), (b) Fixed points: self-consistent solution of K = f{K). 


approach is radically different from traditional statistical mechanical approach 
of constructing ensembles and using probability distribution functions. 

This qualitative picture of Kadanoff ’s was turned into an efficient computa- 
tional procedure with far reaching further physical insights (origin of universality 
etc) by Ken Wilson in early 70s. This approach is known as the renormalization 
group approach of Ken Wilson^]. Although the concept in its primitive form 
and with certain drawbacks was know in the early days of quantum field theory 
(QFT) for example, Gell-Mann and Low etc, it was new approach in critical 
phenomena (in 1970s) that resolved the problems of MFTs in describing critical 
indices. We will clarify the Wilson’s approach using one specific example, the 
</> 4 — theory, although the approach has much more general applicability. The <^ 4 
theory also connects to our previous example of the Ising model! 

Next subsection is devoted to Wilson’s formulation. A brief biographical 
sketch of the man behind this great feat is given in the appendix. 


5.2 Wilson’s formulation 


As promised in the last subsection we would like to answer the following ques- 
tion: How did Wilson change Kadanoff ’s qualitative scheme into a 
powerful computational method? For this we consider the </> 4 - theory. 

The 4 > 4 Theory: Consider again the example of the Ising model (fig- 

ure ED in zero external magnetic field: 


H 

k^T 


K Wj- 

<nn> 


We stress that if one wants to understand the physics of the renormalization 
group one must work through the mathematical apparatus which is rather stan- 
dard now; namely, the functional integrals and the method of contractions. 
Below the logic of the calculation and essential steps are given (see for details 
the original reference; Wilson-Kogut in [3]). The partition function Z for the 
Ising model can be written as: 


p+oo 


/ ds rn 28 {s 2 m -l)e K ^nYji 


S-n Sri+i 


n 


f+OO 


ds r 


*-k bs m- us tn.p K Yl n Yji 


'71 j 
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Smoothing delta functions 





Figure 18: Smoothing the delta function peaks into narrow Gaussian peaks. 


Where on the resulting expression on the right hand side we have smoothed the 
sharp delta function peaks into narrow Gaussian peaks (see figure [18]). Writing 
the above in fourior space with: 


1 


where 


+7 T /»+7T 

dqi.... / dq d , 
J —TT 


transforms the partition function to: 




f+OO 


dsjYiC 


f 'Y' ■ \ e ' qj ~ l| 2 + & — 2/fd)cr q cr_ q — u f f f er q . er q2 Oq-. C-q, - 
j<7 ^ J qi J q2 J Q3 


The program is: We will perform one quick step and then the sequence of 
steps! With this we will end up at two coupled algebraic equations. From there 
we can compute the critical indices. 

One quick step: 

Let us assume that the atomic scale fluctuations do not effect the value of the 
critical indices (an assumption which can only be verified a-posterior). Thus we 
will consider wavelengths much larger than the lattice constant. A momentum 
cut-off |q| < q c will be imposed in the integrals (where q c « -ff, a = 1, a is 
the lattice constant). Thus one can approximate \e lqi — 1| 2 ~ qf. But limiting 
|q| < q c harms: 



One cannot get back the original s rn (as large momenta |g| > q c are excluded in 
the limits of the integral!) So there is a technical problem. To resolve this one 
is forced to define the functional integral to count all possible states that the 
system can realize: 


Z = 


Where 


And now 


2 _/) O + r ) (T q a '-q U J qi J q2 Sq 3 CT< 12 CT< 13 a — <31 — <12 — 33 r ^ 2I\.d 


K 


p+oo 


= n / (•■••)* 

f 1 /‘+ < 3'c f 


Qj 


-\-Qc 


This is one typical example of defining an Effective Field Theory (EFT)! Our aim 
is to execute Kadanoff ’s idea of “thinning out the degrees-of-freedom” . Here we 
are doing in momentum space (instead in real space as was done in the previous 
example of the last subsection). Higher momentum states will be integrated 
out iteratively. The logic of the procedure is given below. The fundamental 
importance of Wilson’s iterative procedure (the renormalization group) is that 
we will get the pivotal function f(K) of the previous example. From that it is 
possible to get the critical point and all the critical indices etc. 


The Procedure: Define the functional integral as: Z = 


with H[a\ = 


I cr 

-h Iq( ( l 2 + r ) a q^-q-w f qi f q2 f q3 cr q2 cr q3 CT_ qi _ q2 _ q3 . This functional integral 

can be written as the product of two functional integrals (as is clear from the 
definition of the functional integral above). 


Z = 


o H[a 0 +(Ji} 


I <ti 


higher momenta integrals 


Here 

_ _ / cro.q = cr 0 */ 0 < |q| < ^ 
q 1 cri.q = cr-l if y < W < 9c 

We want to integrate out oys (the higher momentum shell). Also rename various 
parts of the Hamiltonian as: 


H W\ = ~\ f (O 2 + r )° qcr-q + ~ l 


(T qi lT q2 <T q3 cr -qi“q2-q3 • 


4i •> qi -'43 


Hf[ct\ 


Hi[a\ 


Let the integration of higher momenta states yields: 


e H 'W\ = gffpko] 


e H F [CTi]+ir/[<T 0 +(Tl] 


O’! 


Here we define a’ = a'^ = i(7 0 ,q (be., o' and q' are rescaled variables with 
0 < |g'| < q c and 0 < |g| < ^). Hp[ao] can be written in terms of scaled 
variables trivially as: 


1 

2 


' 0<|q| < - 


(q 2 + r)cr q cr_ 


__ (c 2 2- d - 2 ) 


0<|q'|<4c s 


( q' 2 +4 r) cr'q/cr^q,. 


notice\ 


This trivial scaling cannot be done on the interaction term (due to coupling of 
momenta). One must use perturbation theory and express: 

Where x = — Hj[a\. In this series expansion inside the functional integral one 
will encounter terms like: 


J(qi,q 2 ,....q fc ) = 


cr l,qi cr l,q2 • 


•^Xq/e^ F ^ ^ {^Xqi &1,CL2 ^1^3 ^1^4- 


<7l 


....all combinations...} . 
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Figure 19: (a)Iteration formula for u n . (b) Calculation of v: if one starts with 
some initial ito ; r o at T = T c , after sufficient number of iterations r — > r* , u — > 
u* (the fixed point) . If one starts at T ^ T c iterations take the initial conditions 
to some asymptotic line D. 


For the proof of it see the original article: Wilson-Kogut in [II]. After doing the 
functional integral over <ns (say up to 2nd order) one must rescale the <tq as 
before (by expanding the momenta 0 < |q| < ^ to 0 < |q'| < q c etc). Effective 
Hamiltonian (after integrating upto second order (see for details Wilson-Kogut 
in 0)) takes the form: 


H W] = -\ f W 2 + T-XV q , + -v! f j J < ; <'<^ q ,_ q ,_ q , . 

J a J 9i J q 2 J % 




Hi[a'] 


With 


r' = 4(r + 3c- ) 

v 1 + r J 

u' = 2 4 ~ d (u — 9c- U 


" (1 + r) 2 


)■ 


These are the coupled algebraic equations. One can make analogy to the 
equation = /(/C^/Zj) that we derived by using Kadanoff’s real space 

blocking method. Here calculation is done in momentum space. By invoking 
the same logic that at the critical point the function / becomes scale indepen- 
dent due to the divergence of correlation length. This divergence of correlation 
length leads to fixed point solution (as we learn in Kadanoff’s method in last 
subsection). Thus self-consistent solution of the above contains the physics of 
the critical point. One can obtain the self consistent solution by iterative pro- 
cedure i.e. , by evolving the above equations as a recursive procedure. This 
iterative procedure can be depicted by the following diagram (figure [191 alb 
Let us express the equations after nth iteration as 

u 

f n+i = 4 (r„ + 3c — —) 

1 + r„ 

u n+1 = 2 4 ~ d (u n - 9 c (1+ ^ )2 ). 
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After sufficient number of iterations one will end up at the critical point u*, r* . 
For illustration consider the case of u n (figure UW a) 1 . The two curves (one for 
u and other for f(u)) meet at point u* . If one starts with some initial value uo 
(uq > u* or uo < u*, does not matter) one will end up, after sufficient number 
of iteration, at u*. 

The great advantage of the Wilson’s formulation is that we are able to know 
the required function explicitly. The other important point is that the critical 
fixed point r * , u* does not depend on the initial parameters in the Hamiltonian 
i.e., 7*0, uo- Thus the critical point behavior does not depend on the microscopic 
details of the system. This is called the universality. For example, /3 = 1/3 
both for the magnetic systems and for the liquid gas transitions! (see for details 

000E3 M)- 

This how Wilson changed Kadanoff’s qualitative ideas into a powerful com- 
putational procedure, as we know today, RG has been applied to wide variety 
of problems. Here by considering one example we have covered the basics of the 
method. 

One can also compute the critical indices. Let us take one example of the 
calculation of v. For this we need to linearize the fixed point equations. For T 
very near to T c (say curves 2, or 3 in figure flWbH: 


r / (T)=r z (T c ) + (T-T c )^| To + ... 
u l (T) = u l (T c ) + (T-T c )^\ Tc + ... 

As r = d is a function of temperature through K = p^p, and u = ppp is 
also a function of temperature. Thus these can be Taylor expanded for T near 
T c . For very large l: 


Or 


ri (T)=r* +{T-T c )^\ Tc +... 
ui(T) = u* + (T-T c )^\ Tc + ... 

r i+i{T) — r* = (T — Tc)—^r\ Tc + ■■■ 

Ui+ 1 (T)-u* = (T-T c )^\ Tc + ... 

For T sufficiently close to T c . As r/+i = f(ri,uf) and tq+i = g(ri,ui) 
(our fixed point equations), the temperature derivatives can again be Taylor 
expanded around the fixed point. This leads to: 


n + i - r* = ci(r ; - r*) + c 2 (rtz - u *) 
ui + i -u* = c 3 (r; - r*) + c 4 (zq - iT) 

As the functions f(ri, uf) and g{ri, ui) are known one can work out the expansion 
coefficients explicitly (see for details Wilson-Kogut in ®Y 


rz+i-r* N 

) — M2x2 ( 

ri ~ r * ) 

Ul+i - u 

/ \ 

. Ul-U J 
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By iterating n— times: 


( ri +n - r* 
V u l+n - U* 


( M'2x2) n 


n - r* 
ui - u* 


The purpose of iterating n-times is that the matrix M n is completely dominated 
by the largest eigenvalue A" (A = 4 for u = 0 case). One can find eigensystem 
of the above matrix and finally (as the largest eigenvalue dominates!): 


ri+n - r* ~ A”whq(T - T c ). 


here Ai is the largest eigen value, w\\ is the corresponding eigenvector, and c; is 
a constant (see for details Wilson-Kogut in[3]). And a similar expression holds 
for the ui. One observe that: 

ri +n+ i(T)-r* at (T - T c = -L) = r l+n (T) - r* at T-T c = t 

M 

Thus the correlation length: 

^(n+n+l i ^Z+n+l) | T=T c +r/ Ai • 

Now 

2 ;+1+i ^( T c + t/A 1 ) = £(r c + r ), as ^(n+„, u i+ „) = u 0 ). 

If correlation length scales as: 


£(T C + t) oc r " 


Then from above: 



2r~ 


ln2 

InXi 


If u = 0 (simple Gaussian model), then Ai = 4 which leads to v = 0.5 (the MFT 
result). With explicit calculation of Ai for it ^ 0 (see for details Wilson-Kogut 
in[5]), v comes out to be: 


v ~ 0.58 


This is very close to the experimentally observed value v = 0.6 (in MFT it is 
0.5!!!!). Similarly other critical indices can be calculated and are quite close to 
what is observed experimentally [IT. I27| . 


6 Why do MFTs fail and how does RGT rectify 
that? 

Thus, by exploiting scale invariance at the critical point and by finding the crit- 
ical point from fixed point equations by iterative procedure one is able to calcu- 
lated correct values of critical indices. The important point that goes implicit in 
the above statement is that iterative procedure is not only a mathematical de- 
vice to reach to the fixed point but also at each step of iteration one is averaging 
over smaller length scale (or higher momentum scale) fluctuations and with this 
appropriate averaging over increasingly longer length scales one is able calculate 
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the correct values of critical indices. With this sophistication of renormalization 
group method we immediately see the drawback of MFTs. In MFTs we do an 
abrupt step in dealing with the fluctuations: (a) = tanh(/i + -y-(er)) i.e. , we 
average out only the smaller length scale (of the order of lattice constant) fluc- 
tuations (in calculating the effective field seen by the test spin) . Near the critical 
point just below the critical temperature, ensemble averages over smaller length 
scales (of the order of lattice constant) lead to zero value of the order parameter 
(here mean magnetization). Only when one averages over a large length scale 
(of the order of correlation length) can one see some tiny magnetization just be- 
low criticality. Thus our assumption of tiny mean magnetization over a smaller 
length scale (of the order of lattice constant) in MFT breaks down! And MFT 
predicts wrong values of critical indices ( see also [7]). 

7 Appendix: K. G. Wilson— a short biographical 
sketch 



Figure 20: Kenneth Geddes Wilson (8 June 1936 - 15 June 2013) with Hans 
Betlie and colleagues; 1982 Nobel Prize in Physics for the renormalization group 
method (Image source: http://news.cornell.edu/stories/2013/06/physics-nobel- 
laureate-kenneth-wilson-dies) . 

Ken Wilson earned his Ph.D. from the California Institute of Technology 
(CalTech) in 1961, studying under Murray Gell-Maan. Elementary Particle 
Physics (now called high energy physics) was a fashionable subject in those 
days (today also it is one of the mainstream topics along with condensed matter 
physics) and students would like to opt for it. But Wilson rebelled, and worked 
for a while on plasma physics. Something that did not go very wcll[3]. During 
that time he used to discuss his problems with S. Chandrasekhar (cousin of sir 
C. V. Raman). Here is an interesting story [T21: 

His father had told him: When you go to Caltech, be sure to meet the two 
great physicists there namely, Feynman and Gell-Mann’. K. Wilson says that 
when he called on Feynman, the latter seemed to be gazing at the ceiling. He 
asked Feynman, ’Professor, what are you working on at present?’ Feynman 
replied, ’Nothing! Wilson then went to Gell-Mann and became his student!!! 

i / 2 \ 

He worked on the solutions of Gell-Mann and Low equation A 2 ^ = 

(f>{e \ , m 2 /X 2 ) and on fixed-source meson theory for his PhD Thesis. Around 
1963 he read Onsager’s paper and became interested in critical phenomena, 
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and after that studying other works on critical phenomena he developed the 
renormalization group theory, using tricks and techniques that he had learnt in 
particle physics. 
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